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Abstract Space plasmas are collisionpoor and kinetic effects prevail leading to wave 
fluctuations, which transfer the energy to small scales: wave-particle interactions 
replace collisions and enhance dispersive effects heating particles and producing 
suprathermal populations observed at any heliospheric distance in the solar wind. 
At large distances collisions are not efficient, and the selfgenerated instabilities con- 
strain the solar wind anisotropy including the thermal core and the suprathermal 
components. The generalized power-laws of Kappa-type are the best fitting model 
for the observed distributions of particles, and a convenient mathematical tool for 
modeling their dynamics. But the anisotropic Kappa models are not correlated with 
the observations leading, in general, to inconsistent effects. This review work aims 
to reconcile some of the existing Kappa models with the observations. 



1 Introduction 

Direct in-situ measurements in the solar wind and terrestrial magnetosphere indi- 
cate that the velocity distribution functions (VDF) of space plasma particles are 
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quasi-Maxwellian up to the mean thermal velocities (the core component), while 
they exhibit non-Maxwellian suprathermal tails (the halo component) at higher en- 
ergies (see recent reviews of Marsch (2006) [1] and Pierrard & Lazar (2010)(2|, and 
references therein). Processes by which the suprathermal particles are produced and 
accelerated |3]|4][5][6][7]|8) are of increasing interest for applications in astrophysics 
and laboratory or fusion plasma devices where they are known as the runaway par- 
ticles decoupled from the thermal state of motion. The solar wind generally appears 
to involve an abundance of suprathermal electrons and ions observed to occur in the 
interplanetary medium, and their analysis provides valuable information about their 
source, whether it is in the Sun or outer heliosphere. 

Accelerated particles (including electrons, protons and minor ions) are detected 
at any heliospheric distance in the quiet wind as well as in the solar energetic particle 
(SEP) events associated to flares and coronal mass ejections (CMEs) during solar 
maximum (see reviews by Lin (1998)[9| and Pierrard & Lazar (2010)[2|). A steady- 
state suprathermal ion population is observed throughout the inner heliosphere with 
a VDF close to ~ v~ 5 ifTOl . and, on the largest scales, the relativistic cosmic -ray gas 
also plays such a dynamical role through the galaxy and its halo ifTTI . 



2 Wave instability constraints of the suprathermal anisotropy 

Nonthermal features and kinetic anisotropics like temperature anisotropies, heat 
fluxes or particle streams are also a characteristic of the solar wind and near the 
Earth's magnetosphere UJ. At low altitudes in the solar wind, velocity distributions 
of plasma particles regularly show an excess of kinetic energy transverse to the lo- 
cal mean magnetic field (7j_ > 7ji, where _L, and || denote directions with respect 
to the magnetic field) most probably due to the compression exerted by the strong 
guiding magnetic field near the Sun. At larger heliospheric distances, the anisotropy 
is controlled by the Chew-Goldberger-Low mechanism: the adiabatic expansion of 
the poorcollision plasma increases the pressure and temperature along the magnetic 
field leading to 7ji > T±. In more violent interplanetary shocks resulting after solar 
flares and coronal mass ejections (CME), injection of particle beams into the ionized 
interplanetary medium creates additional anisotropy [Q]|9]. 

The high rate of occurrence of an excess of perpendicular temperature (T± > Tj|) 
in measurements at large distances lfT3] [141 is a proof that other mechanisms of ac- 
celeration, namely, the wave-particle interaction, must be at work there dominating 
the adiabatic expansion. Indeed, large deviations from isotropy quickly relax by the 
resulting wave instabilities, which act either to scatter particles back to isotropy, or 
to accelerate lower energetic particles (Landau or cyclotron heating) and maintain a 
suprathermal abundance because thermalization is not efficient at these scales JTJ. 

None of these processes is well understood, mostly because these plasmas are 
low-collisional and require progress in modeling the wave turbulence, going beyond 
MHD models to use a kinetic and selfconsistent description. In such plasmas trans- 
port of matter and energy is governed by the selfcorrelation between particles and 
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Fig. 1 Electron velocity dis- 
tributions observed by Wind 
mission at 1 AU, as energy 
spectra (top) parallel (solid) 
and perpendicular (dashed) to 
magnetic field; and velocity 
space contours (bottom) in a 
high-speed solar wind. Note 
the anisotropic isodensity 
contours, less for the core 
and more pronounced for the 
halo and the strahl in fast 
solar wind (after Pierrard et 
al. (1999) fT2l ). 
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electromagnetic fields, which can, for instance, convect charged particles in phase 
space but are themselves created by these particles. Thus, the suprathermal popula- 
tions involve selfconsistently in both processes of wave turbulence generation and 
particle energization, and the resulting Kappa functions that elegantly describes dis- 
tributions measured in the solar wind J2], represent therefore not only a convenient 
mathematical tool, but a natural and quite general state of the plasma lfl5l . 
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Fig. 2 Solar wind electron temperature anisotropy (7j_/7j|) versus J3||: the instability thresholds 
from a bi-Maxwellian model do not constrain the suprathermal halo (after Stverak et al. (2008)). 



While the thermal core in the solar wind is less anisotropic, the superthermal 
halo has in general a pronounced temperature anisotropy (7j_ ^ 7j|) with respect 
to the local magnetic field 0][2]- Deformations of the VDFs observed in the so- 
lar wind are not as strong as one would expect from a free motion of particles 
(Kasper et al. 2002, Stverak et all. 2008, Bale et al. 2009). Because collisions are 
not effective, any increase of kinetic anisotropy is limited by converting the excess 
of free energy into electromagnetic fluctuations, and pitch-angle diffusion. Model- 
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ing the VDF with a bi-Maxwellian, the instability thresholds shape precisely the 
core (Fig. 2, left) but show regularly an important departure from to the halo limits 
(Fig. 2, right). While the whistler instability limits the perpendicular temperature 
to grow (Tjl > Tii), the firehose instability constrains any excess of parallel tem- 
perature (7j_ < 7j|). Suprathermal particles cannot fit into a Maxwellian approach 
( K — > °°), but must be modeled with an appropriate Kappa function. A bi-Kappa (or 
bi-Lorentzian) function has extensively been used to model gyrotropic distributions 
and their dispersion/stability properties [jT6] [T7] 

i r[K+i] ( v j v 2 , \~ K l 

Fx V||,vj_ =^72~2 TT^T 1 + ^ + ^ ■ CD 




Fig. 3 Firehose instability: thresholds (left) and growth rates (right) from a bi-Kappa model. 
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Fig. 4 Whistler instability: thresholds (left) and growth rates (right) from three different models: 
the bi-Kappa, the product-bi-Kappa and the bi-Maxwellian (after Lazar et al. (201 la)). 



However, this model does not provide the expected better fits to the observations. 
For instance, both instabilities of interest for the electron populations, the whistler 
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and the firehose, are inhibited and threshold constraints do not approach but rather 
depart from the measured limits of the halo |[T8] [191 |2T1 1221 . Marginal conditions of 
the firehose instability are illustrated in Fig. 3 (left) showing the need for a larger 
temperature anisotropy and a larger plasma J3|| = SnnksT^ j B\ to produce the in- 
stability in Kappa distributed plasmas (low values of fc — > 3/2). By comparison to 
a bi-Maxwellian (fc — s- °°), the growth rates are in general reduced and restrained 
to small wavenumbers (Fig. 3, right) |[T9l l22l . The same tendency is observed in 
the evolution of the whistler instability (Fig. |4j: by comparison to a bi-Maxwellian 
thresholds do not depart much, but the growth rates decrease significantly for lower 
values of the index K (Fig.g] right) lfl"8ll2D . 



3 Product-bi-Kappa model 




It is evident that a novel Kappa function is needed to model kinetic anisotropics and 
to incorporate the excess of free energy expected to exist in suprathermal anisotropic 
plasmas providing better fit to the observations. Recently a new approach has been 
proposed 11201 1211 |23l based on the anisotropic product-bi-Kappa function |[T6l 

F 2(n^) = ^itr ,/;_" ' J „ 1 1 + ^75 1 " + - (2) 

By contrast to the bi-Kappa function |T|), which seems to be less realistic because 
the two degrees of freedom are coupled and controlled by the same power index 
fc, the new product-bi-Kappa function shows an advanced flexibility in modeling 
gyrotropic VDFs with two distinct temperatures 7ju and two distinct power indices 
Kji j_. Thus, a new concept for the particle anisotropy can be introduced by including 
both anisotropies of the temperatures 7j| ^ T± and the Kappa indices Kj| ^ K^_. The 
analysis becomes therefore more relevant but complicated and this is probably the 
reason this model was only occasionally invoked (after Summers & Thorne (1991) 
proposed it), merely in a simplified Maxwellian-Kappa form ifTTl 

«(,,,„) = ^'j f , + X] (-4) • 

11 7Z i ' 2 W\\ K l ' 2 r[K+ 1/2] y KW \\J V W ±J 

The distribution function (01 reduces to distribution function (O in the limit of very 
large values of K± — > °°. Thus both forms of the new model represented by the 
general distribution function (|2) or by the particular form (0 reduce to the same 
bi-Maxwellian in the limit of very large power indices (Kni — > °°)- 

Contour plots in velocity plane are illustrated in Fig. q\ showing no visible ex- 
cess of asymmetry of the bi-Kappa (dotted-dashed lines) by comparison to the bi- 
Maxwellian (solid lines), but a prominent asymmetry and anisotropy of the new 
product-bi-Kappa distribution function (dashed lines) by comparison to both the bi- 
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Fig. 5 Two sets of contours, 
one close to maximum and 
one close to the minimum 
base: the bi-Maxwellians 

= K± -> oo) with solid 
lines, product-bi-Kappas 
with red dashed lines, and 
bi-Kappas with blue dotted- 
dashed lines (for the same 
Vj i| /c = 2vr t x/ c = 0.2, and 

K|| = K± = 3). 



V\\/C 



Maxwellian and bi-Kappa distribution functions, even for the same temperatures 
7] I = T±, and the same Kj| = K± = K ll20ll2Tl . However, the new distribution model 
and its dispersion properties and stability must directly be confronted to the obser- 
vations in the solar wind and terrestrial magnetosphere. 





Fig. 6 The electron VDF with a strahl component in the high speed solar wind at a distance (in 
AU) R =0.98 (left), 0.64 < R < 0.70 (midle) and 0.29 <R< 0.34 (right) (after Pilipp et al. 1987a). 



Thus, contours of the electron distribution functions measured at different helio- 
spheric distances in a high speed solar wind (see Figs.|6]and|7J show such skewed, 
highly anisotropic tails along the magnetic field direction ||27ll28l , which look quite 
similar to the product-bi-Kappa contours in Fig. In Fig. [6] skews of the electron 
distributions are asymmetric being produced by the magnetic field aligned strahl 
population, which is highly energetic (suprathermal) and antisunward moving |[27l . 
The presumable origin of the strahl electrons is in the energetic ejecta from the coro- 
nal holes, and radial evolution show a decreasing with distance from the Sun in the 
favor of the halo population, which is enhancing 1291 . The strahl component is the 
main driver of the (electron) heat flux, and the main contributor to the anisotropy of 
suprathermals, apparently a manifestation of the adiabatic focusing. 
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The quite-time distribution in the fast solar wind (Fig. |6j is modeled by the 
product-bi-Kappa function only on one side, namely, the outward direction (vii > 
0), while semicontours in the backward direction are similar to a bi-Kappa (low 
anisotropic) model. Analytically one can combine the Eqs. (O and (01 to form 



F 4 (v||,vl) = H[-v.|]C L 1 
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where H[x] is the Heaviside function, and Cl.r are normalization constants. Contour 
plots are similar to the one-sided strahl observations in Fig. [6] Hence the immediate 
effect would be a corresponding asymmetric spectrum of the selfgenerated wave 
fluctuations, which are expected to dominate the outward direction. The future work 
should explore the stability of this new model and correlate with the observations. 

Electron distribution functions (see Fig. |7J with an unusual double strahl have 
been observed on rare occasions in the solar wind by Helios probes |f28l . Contours 
plots of these distribution function are illustrated in Fig . [7] (left) and show a strong 
symmetric bidirectional anisotropy very similar to the product-bi-Kappa model in 
Fig. [5] Double-strahl distributions would be observed if the spacecraft was fortu- 
itously near the outer end of a magnetic loop connected to the Sun, and electrons 
coming from both foot points have traveled roughly the same distance and thus may 
form a symmetric distribution with a counterstreaming strahl aspect. A similar dis- 
tribution is expected in the observations if the spacecraft was at the outer end of a 
disconnected large scale loop ll28l . 




Lug IF/f„„l 
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Fig. 7 Contours and one-dimensional cuts of unusual double-strahl distribution function observed 
by Helios (after Pilipp et al. 1987b). 



Another important confirmation should come from the resulting electromagnetic 
fluctuations which constrain particle velocity anisotropy to grow |[T4l l24l . Some 
preliminary analysis have indeed shown that the whistler instability growth rates 
(see Fig. [4] right) calculated for the general product-bi-Kappa model are enhanced 
(dashed line) by comparison to those obtained for a bi-Maxwellian (solid line) or a 
bi-Kappa model (dotted line), and the instability thresholds illustrated in Fig.|4](left) 
tend to approach better the limits of the electron halo in Fig . [2] (right) ll2D . Unlike 
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the bi-Maxwellian or bi-Kappa distributions, the product-bi-Kappa anisotropy can 
be stable against the excitation of the low-wavenumber whistler waves, the critical 
wavenumber£ c . ~ £1/{B\\ */Kj|A) is vanishing only for K\\ — > °° (Maxwellian plasma). 
Moreover, it was shown the instability of an asymmetric Maxwellian-Kappa distri- 
bution without an effective temperature anisotropy (7j_ = 7j|) against the excitation 
of the whistler waves propagating obliquely to the magnetic field 11301 , most proba- 
bly due to the same asymmetry of contours in velocity space (Fig.[5]l. Thus, the new 
product-bi-Kappa model seems to incorporate the excess of free energy expected to 
exist in anisotropic suprathermal plasmas. For the opposite case ( 7j| > T±), the fire- 
hose instability is currently under active investigations, and, according to the preced- 
ing discussion, we expect that a product-bi-Kappa model will enhance the growth 
rates and provide better agreement of the marginal stability with the observed limits 
of the halo (Fig. [2] right). 

Generation of such asymmetric distribution profiles is also supported by the 
anisotropic turbulence of the solar wind, where the acceleration is expected to occur 
either perpendicular (cyclotron damping)|]5] or in direction (Landau damping)|6| 
of the interplanetary magnetic field. The nonlinear wave-wave and wave-particle 
couplings involving intense low-frequency Alfven waves or electrostatic Langmuir 
and ion sound (weak) turbulence driven by the beam-plasma instabilities can also 
be responsible for the acceleration in the corona, solar wind and magnetosphere 
|01|7l. However, these models do not provide much guidance on a full 3D evolution 
of the distribution function and its nonthermal features, like suprathermal tails and 
kinetic anisotropics, in the acceleration process. This task is still complicated, but 
will, hopefully, make the object of the next investigations. 



4 Conclusions and perspectives 

A large variety of nonthermal features like temperature anisotropics, heat fluxes or 
particle streams are permanently observed in the solar wind and near the Earth's 
magnetosphere 0] [2]. These do not grow indefinitely, but tend to regulate them- 
selves since they contain sufficient free energy to drive plasma instabilities and 
micro-turbulence 11251 . Plasma wave turbulence provides the main dissipation mech- 
anisms maintaining a fluid-like behavior of the plasma in spite of the fact that col- 
lisional free paths become comparable to the large scales of the heliosphere. Bi- 
Maxwellian models have extensively been used to describe particle distributions 
and their dynamics in the solar wind, but the observations clearly indicate the rele- 
vance of Kappa distribution functions which incorporate both the quasithermal core 
and the suprathermal halo. Because the anisotropic Kappa functions are not corre- 
lated, in general, with the observations, here we have reviewed these models, and 
established their direct or indirect confirmations by the observations. 

The nondrifting distribution functions used to describe the temperature anisotropy 
of the different plasma components (with respect to the magnetic field), are the bi- 
Kappa and product-bi-Kappa functions. The bi-Kappa function seems less realistic 
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because the two degrees of freedom parallel and perpendicular to the magnetic field 
are coupled and controlled by the same power index k. Instead, the new product- 
bi-Kappa function shows more flexibility in modeling the gyrotropic VDFs with 
two distinct temperatures T»j_ and two distinct power indices K\\ t j_. A simple com- 
parison of the contours plots show indeed a significant excess of anisotropy of the 
product-bi-Kappa by comparison to both the bi-Kappa and bi-Maxwellian models. 

A further comparison with the contours plots of the measured distributions indi- 
cates two important results of our analysis. First, the bi-Kappa model is appropri- 
ate to describe particle distributions in the quiet and slow solar wind with a min- 
imum strahl influence and, in general, small deformations of the halo. Secondly, 
the product-bi-Kappa (including the Maxwellian-Kappa) function seems to be ad- 
equate for modeling particle distributions in the fast solar wind with a prominent 
strahl component in the direction of magnetic field. Moreover, a double-strahl dis- 
tribution looking just like the product-bi-Kappa contours can be observed near the 
magnetic large scale loops with electrons counterstreaming from both foot points. 

The selfgenerated wave spectra originating from the Kappa models exhibit dif- 
ferent properties. Thus, the wave instabilities constrains calculated for a bi-Kappa 
model do not fit with the halo limits in a slow solar wind as was expected, but those 
resulting from a product-bi-Kappa model seem to shape the halo electrons for any 
slow or fast solar wind. This suggests that the new product-bi-Kappa model incor- 
porates in a proper way the excess of free energy expected to exist in anisotropic 
suprathermal plasmas. Assuming that the same wave fluctuations are the most plau- 
sible mechanism of acceleration and formation of suprathermal populations, we 
conclude pointing out the relevance of our study by the fact that an appropriate 
Kappa model enables a powerful selfconsistent analysis as it should be itself a prod- 
uct of particle acceleration by the selfexcited fluctuations. 
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